Relativistic Hartree-Bogoliubov theory in coordinate space: finite element solution for a nuclear system with spherical symmetry

A C++ code for the solution of the relativistic Hartree-Bogoliubov theory in coordinatespace is presented. The theory describes a nucleus as a relativistic system of baryons andmesons. The RHB model is applied in the self-consistent mean-field approximation to thedescription of ground state properties of spherical nuclei. Finite range interactions areincluded to describe pairing correlations and the coupling to particle continuum states.Finite element methods are used in the coordinate space discretization of the coupledsystem of Dirac-Hartree-Bogoliubov integro-differential eigenvalue equations, and Klein-Gordon equations for the meson fields. The bisection method is used in the solution of theresulting generalized algebraic eigenvalue problem, and the biconjugate gradient methodfor the systems of linear and nonlinear algebraic equations, respectively. Work supported by GSI Darmstadt (GSI TM DIT 36-46-4030)E-mail:[email protected] von Humboldt Fellow, on leave of absence from University of Zagreb, Croatia1